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ABSTRACT 

Context. The solar photospheric abundance of oxygen is still a matter of debate. For about ten years some determinations have 
favoured a low oxygen abundance which is at variance with the value inferred by helioseismology. Among the oxygen abundance 
indicators, the forbidden line at 630nm has often been considered the most reliable even though it is blended with a Nil line. In 
Papers I and Paper II of this series we reported a discrepancy in the oxygen abundance derived from the 630 nm and the subordinate 
[O I] line at 636 nm in dwarf stars, including the Sun. 

Aims. Here we analyse several, in part new, solar observations of the the centre-to-limb variation of the spectral region including the 
blend at 630 nm in order to separate the individual contributions of oxygen and nickel. 

Methods. We analyse intensity spectra observed at different limb angles in comparison with line formation computations performed 
on a C05BOLD 3D hydrodynamical simulation of the solar atmosphere. 

Results. The oxygen abundances obtained from the forbidden line at different limb angles are inconsistent if the commonly adopted 
nickel abundance of 6.25 is assumed in our local thermodynamic equilibrium computations. With a slightly lower nickel abundance, 
A(Ni) ~ 6.1, we obtain consistent fits indicating an oxygen abundance of A(0)=8.73 ± 0.05. At this value the discrepancy with the 
subordinate oxygen line remains. 

Conclusions. The derived value of the oxygen abundance supports the notion of a rather low oxygen abundance in the solar photo¬ 
sphere. However, it is disconcerting that the forbidden oxygen lines at 630 and 636 nm give noticeably different results, and that the 
nickel abundance derived here from the 630 nm blend is lower than expected from other nickel lines. 

Key words. Sun: abundances - Stars: abundances - Hydrodynamics 


1. Introduction 

Although oxygen is the third most abundant element in the 
Sun (after hydrogen and helium), its pre cise abund a nce is cur¬ 
rently under discussion. Very early on, ISedlmavJ (Il974t) de¬ 
rived a rather low oxygen abundance in comparison to the 
accepted value at the time, and attributed the difference to 
the important role of departures from local thermodynamic 
equilibrium (LTE) for the oxygen infrared triplet at 777 nm. 
Th e debate gained mome n tum th a nks to the pioneerin g work s 
of iKiselrnan & Nordlund (1995h. Allende P rieto et al.l (l200lh . 
iHolwegerl (IZOOlll . and Asnlund et al.l (l2004l) who were the first 
to apply hydrodynamical simulations (hereafter 3D models) to 
the determination of the solar oxygen abundance. Their results 
advocated a downward revision of the solar oxygen abundance, 
and therefore of the solar metallicity Z, whi ch now appears to be 
in conflict with the helioseismic data (see iBasu & Antiall2008l 
for a review of the topic). In the last ten years much effort 


has been dedicated to confirming or dism issing the downward 
revisi o n of the sola r oxygen abundance (lAsplund et alJ 12004 
I2OO9I1 . iPereira et alJ (l2009bh analysed the [OI] line; they put 
A(Ni) =6.22 and derived A(0)-8.64 from their centre-disc spec¬ 
trum. Ayreset^ J2006ll investigated the infrared lines of CO 
while Avres I (l2008h analysed the [O il line at 630 nm, both in- 
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dicatin g a “high” oxygen abundance. ICenteno & Socas-Navarr^ 
(l2008h also analysed the [OI] line at 630 nm and derived a high 
value, A(0) = 8.86+0.07 ]Socas-Navarro & Norton|(l2007h anal- 
ysed t he oxygen triplet deri ving a “low” NLTE oxygen abun¬ 
dance. ISocas-Navarrol (12015h investigated the [OI] line with the 
support of an empirical three dimensional model, obtaining a 
high A(0), but with a large uncertainty and with the possibil¬ 
ity of deriving a low oxygen abundance by slightly changing the 
model. This paper is particularly interesting because it shows the 
limiting factors of the models we use. 

Our project aims at investigating the oxygen indicators in the 
solar spectrum and deriving the solar oxygen abundance based 
on a CO^BOLD 3D mod el of the solar atm osphere. In a previ¬ 
ous paper of this series (ICaffau et al.ll2008l) . hereafter Paper I, 
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we derived a measure of the solar oxygen abundance from an 
ensemble of atomic lines. Among the lines employed in that in¬ 
vestigation the two [OI] lines at d = 630 nm and 636 nm were 
singled out as potentially problematic. The abundance we de¬ 
rived from the two lines differed by 0.1 dex, both when using the 
3D model and when using the ID Holwege r-Miiller solar model 
(lHolwegeill967tlHolweger & Mulleill974 hereafter HM). This 
is certainly a good reason to investigate these lines further. In 
addition, the [OI] lines, and in particular the stronger 630 nm 
line, have often been considered the best indicators of the oxy¬ 
gen abundance, given that they are weak and formed in condi¬ 
tions of LTE. When using other atomic lines to determine the 
solar oxygen abundance, one is forced to introduce an additional 
uncertainty related to the computation of the NLTE corrections, 
which stems from our lack of knowledge of the efficiency of 
excitation and ionisation of oxygen by collisions with hydrogen 
atoms. Only for a few atoms do we know the cross sections of in¬ 
elastic collisions with neutral hydrogen from either experiments 
or quantum-m echanical computations. The use of the Dra win 
approximation (lDrawinll969 ; Steenbock & Holwegerll98^ has 
been seriously criticised by Barklem et al. ( ffllll). E or further 
discussion on this point see also Steffen et al.l (l20L§) . We esti¬ 
mated the abundance uncertainty due to this unknown parameter 
as about 0.03 dex (Paper I). If one could reliably determine the 
oxygen abundance from the forbidden lines, the role of hydro¬ 
gen collisions could be constrained by requiring that forbidden 
and permitted lines yield the same abundance. 

We investigated in ICaffau et al.l (1201 3h (Paper II) the forbid¬ 
den lines in a sample of stars from high-quality spectra and con¬ 
cluded that the disagreement in the oxygen abundance derived 
from the two [OI] lines is generally evident in dwarf stars, while 
it vanishes in giants. 

In the present investigation we exploit observations of solar 
intensity spectra at different limb angles to put additional con¬ 
straints on the abundance analysis of the forbidden oxygen line 
at 630 nm. 


2. Atomic data and model atmospheres 

Eor the forbidden oxygen line and for the Ni i line we employ the 
same atomic data as used in Pa per I. The adopted /-values are 
from IStorev & ZeippenI (l2000h . The damping parameters have 
no influence on the profile of these weak lines. 

Our analysis is based on a 3D model atmosp here computed 
with the CO^BOLD code (lErevtag et al.l l20l'^ . The box-size 
of the CO^BOLD 3D solar model atmosphere is 5.6 x 5.6 x 
2.25 Mm, with a resolution of 140 x 140 x 150. The simula¬ 
tion covers a range in time of 2.1 h, represented by 20 snapshots; 
this time interval corresponds to 25 periods of the solar 5 min 
oscillations and about 18 convective turn-overtimes. The snap¬ 
shots were taken at equal time intervals and were selected to 
provide a good representation of the statistical properties of the 
complete temporal sequence o btained in the CO^B OLD run. The 
model has been described in ICaffau et al.l (1201 ih and used for 
the chemical investigation of the photospheric solar abundance 
of several elements. The spectrum synthesis code employed is 
Linfor3E[3. The spectra are calculated for the same ji inclina¬ 
tions as the observed spectra, and for four values of azimuthal 
angles. The spectrum corresponding to each /r-value is an aver¬ 
age over all surface elements, all azimuthal directions, and all 
time snapshots. 


' http://www.aip.de/Members/mstefFen/linfor3d/linfor_3D_manual.pdf 


3. Observations 

We analysed five sets observations from 

- Hinode (five disc positions); 

- the Swedish Solar Telescope, SST (four disc positions); 

- the McMath-Pierce Solar Telescope, Kitt Peak (six disc po¬ 
sitions); 

- the Vacuum Tower Telescope, VTT (five disc positions); 

- THEMIS (eight disc positions). 

We think that all of these observations are extremely useful to 
the discussion of the solar oxygen problem. However, two sets 
clearly stand out for different reasons: the Hinode observations, 
taken from space, are not contaminated by any telluric feature; 
the VTT data are of higher spectral resolution than any other 
data set. The Kitt Peak and the SST observations are of simi¬ 
lar quality, although the latter ones have higher signal-to-noise 
ratios. The THEMIS observations present some difficulty in the 
determination of the continuum. In the following, each set of 
observations is discussed in a dedicated section. 

We restrict our analysis to the spectra at // > 0.2 for all 
observed sets of spectra for the following reasons. At smaller 
ju-values sphericity effects become large enough to influence 
the computed line shapes. The most extreme limb spectra are 
formed in more external regions of the solar photosphere, and 
more so for oxygen than for nickel (see Eigs. [T]and|2l respec¬ 
tively). The thermal structure and velocity fields of these external 
low-density layers of the atmosphere are more difficult to model 
correctly, since they are more sensitive to the upper boundary 
condition of the simulation box. 


3.1. Hinode data 

Hinode (JAXA/NASA) observations were obtained with the 
Solar Optical Telescope (SOT, 50 cm aperture) and the Spectro 
Polarimeter (SP) using a thin slit of 0.3 arcsec x 160arcsec. 
Details of the solar granulation are clearly visible (granules, in¬ 
tergranules). Observations were performed at ten locations along 
the solar equator with the slit orthogonal to the equator direction 
so that pixels along the slit are seen at a mostly constant helio¬ 
centric direction, cosine ju. In fact, there is a small fj. variation 
along the slit because pixels at the end of the slit (+80") do not 
exactly have the /i of the slit centre, but it is a small effect. The 
resolution provided by the spectrograph is about 200 000 so that 
the spectra exhibit resolved features of 3 pm at 630 nm. The sam¬ 
pling is 2 pm. The SP is usually centred on the Eei 630.15 and 
630.25 nm lines, but it is possible to move the spectral coverage 
slightly to [OI] 630.03nm and Eei 630.15nm. The CCD pro¬ 
vided spectral images of 112 x 512 pixels corresponding to a 
spectral/spatial field of 0.23 nm x 160 arcsec. Individual images 
have a signal-to-noise (S/N) ratio of 1000 in the continuum. Eor 
each position along the equator, a small scan along 24 arcsec was 
made in the equatorial direction, which is the reason why we se¬ 
lected only the central part of the scan to stay at constant // along 
the equator. Two sides of the CCD were exposed (1.6 s exposure 
time) and we used both. All Stokes parameters were recorded 
(full Stokes mode I, Q, U, V), but we selected only intensities, 
I. However, the full Stokes mode allowed us to check on the Ee i 
630.15nm line (Lande factor 1.67) that only a few pixels were 
affected by magnetic fields of the solar network. The curvature 
of lines was corrected by a parabolic fit derived from the strong 
Eei 630.15nm line close to the [OI] 630.03nm line with flat- 
held data. Then, a summation was done in the spatial direction 
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along the slit to get high precision spectra for various //-values 
with a final nominal S/N ratio of 20 000 in the continuum. 

3.2. SST data 

These spectra were observed at the Swedish Solar Te lescope in 
May 2 007. They are described in iPereira et al.l (l2009ah and anal¬ 
ysed in iPereira et al.l (l2009bl) . The observations were taken with 
slit width of 0.11 arcsec, the spectral resolving power is 200 000, 
the four //-values used are 0.999, 0.793, 0.608, and 0.424. 

3.3. McMath-Pierce data (WCL) 

The observational data consists of solar inten¬ 
sity spectra for nine heliocentric angles 6 (// = 

1.00, 0.87, 0.66, 0.48, 0.35, 0.25, 0.17, 0.12, 0.10), observed 
with the McMath-Pierce Solar Telescope at Kitt Peak, using the 
rapid scan double-pass spectrometer (see lLivingston et ^1120071 
and references therein) in single-pass mode. The observations 
were taken on September 14, 2006; the local time was about 
09h30m to 10h30m. The resolution was set by the slit width: 
0.06 mm entrance, 0.1 mm exit. This yields a spectral resolving 
power of about 90 000. The scan length, in wavelength, is 2048 
points. The image scale is 275 mm '. The slit was set parallel to 
the nearest limb. 

The setting was fixed in order to observe the two oxygen 
forbidden lines at 630 nm and 636 nm. The observations were 
performed on the north-south axis in the sky, from solar centre 
to limb. For each line we obtained two spectra for any //-value 
off-centre, and four spectra of the intensity at disc-centre. Each 
observation freezes the solar five-minute oscillations, covering 
a shorter interval of time. Effective exposure time was 100 s for 
// > 0.25 (ten scans of 10 s), and 200 s for // < 0.17 (twenty 
scans of 10s). We summed the spectra at equal // to improve 
the signal-to-noise ratio. No wavelength calibration source was 
observed and we used the available telluric lines to perform the 
wavelength calibration. The positions of the telluric lines were 
measured on the high-resolution, hi gh signal-to-noise ratio , so¬ 
lar spectra of the disc-centre atlas of Neckel & Labsl (Il984l) . as¬ 
suming a linear dispersion relation. We will refer to this set of 
observations as WCL data. 

3.4. VTT data 

The VTT data were obtained on Tenerife, Canary Islands, 
with the new prototype instrument LARS (Lars is an Absolute 
Reference Spectrograph). The set-up is based on the high- 
resolution vertical Echelle spectrograph of the VTT and employs 
a laser frequency comb (EEC) wavelength calibration source 
(see e .g. lMumhv et al.l2007ilWilken et al.l20r^lSteinmetz et all 
l2008h and a single-mode fibre feed (SMEE). The wavelength cal¬ 
ibration obtained from the EEC is accurate to 1 m s '. Light from 
the VTT focal plane is coupled to the fibre from a circular field 
of view (EOV) with a diameter of 10 arcsec on the sky. Elat fields 
are obtained with a fibre-coupled continuum lamp, resulting in 
an excellent gain-table correction. The set-up provides a spectral 
purity and fidelity that is, to the best of our knowledge, not at¬ 
tainable with any other instrument currently available for solar 
observations. An ea rly versio n of the instrument is described by 
iDoerr et ^ (I2012alf9) and bv iProbst et al.l (l2015h . 

The spectra used in this work were obtained during several 
campaigns in the testing and commissioning phase of the instru¬ 
ment in August and December 2013 and in March and July 2014. 


Table 1. Census of the VTT data. Data from several individual 
observations is summed for each heliocentric angle, correspond¬ 
ing to total observation times between 36 and 228 minutes. 



Nobs 

Tobs (min) 

S/N 

1.0 

8 

228.4 

6900 

0.8 

2 

49.9 

4600 

0.6 

5 

127.6 

6400 

0.4 

2 

35.8 

3200 

0.3 

2 

53.6 

4400 


Repeated observations were carried out at disc-centre and at he¬ 
liocentric angles corresponding to // = {0.8,0.6,0.4,0.3,0.2}, 
but only the data for // = 1.0 .. .0.3 were used. Individual ob¬ 
servations ran between 15 and 30 minutes, depending on the ob¬ 
serving conditions. Seeing usually was moderate to bad which is, 
however, not important as the instrument was fed with the spa¬ 
tial average of a 10 arcsec EOV in all cases. We took some care 
to avoid regions with magnetic activity by monitoring G-Band 
context images and Calcium II K full-disc images provid ed by 
the ChroTel telescope in real-time (iKentischer et al.ll2008h . The 
absolute pointing accuracy of the VTT is better than 5 arcsec and 
was checked frequently and adjusted if necessary. 

All spectra were corrected for dark-current, and the gain- 
table obtained from the continuum lamp was applied. Residual 
spectral slopes were corrected by fitting a second-order polyno¬ 
mial to selected high-points in our spectra and to the correspond¬ 
ing points in the Kitt Peak Eourier Tran sform Spectro meter 
(ETS) disc-centre spectrum as provided bv lNeck^ (Il999li . The 
ratio of both polynomials gave the correction factor that was ap¬ 
plied to our spectra to force their continuum to be identical with 
that of the ETS spectrum. The wavelength calibration was de¬ 
rived from the EEC calibration spectra. The laser lines are com¬ 
pletely unresolved by the spectrograph and can be used to mea¬ 
sure its instrumental profile (IP). Prom the full width at half max¬ 
imum of the IP, we derived a spectral resolution between 600 000 
and 750000 for these observations, depending somewhat on the 
alignment of the SMEE, which was changed between the indi¬ 
vidual observing campaigns. 

The calibrated spectra were interpolated to an equidistant 
wavelength scale of ImApixeE' (native sampling is 2.73 mA 
pixel '). The EEC calibration naturally removes the effect of 
any instrumental drift and the radial velocity of each observed 
spectrum was subtracted, resulting in a wavelength scale that is 
defined in the solar reference frame. All spectra measured at the 
same heliocentric angle were then summed to increase the S/N 
ratio. The S/N ratio was measured in the continuum at 6300 A 
and it reaches a value of 6900 for the disc-centre spectra. Table[T] 
provides a census of the calibrated spectra that are used in this 
work. 

3.5. THEMIS data 

THEMIS observations were obtained with the eight-metre spec¬ 
trograph in classical spectroscopic mode using a thin slit of 
0.5 arcsec x 120 arcsec. However, the seeing was never better 
than 1 to 2 arcsec, so that details of the solar granulation are 
not visible. Several scans along the solar radius were performed 
with the slit orthogonal to the radius direction (so that pixels 
along the slit are seen at almost constant //). The interference or¬ 
der was 36 and was selected by the predisperser spectrograph. 
The final resolution provided by the eight-metre spectrograph 
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is about 300 000 so that the spectra exhibit resolved features of 
2 pm at 630 nm. The CCD provides spectral images of 512x512 
pixels corresponding to a spectral/spatial field of 0.5 nm x 120 
arcsec with exposure time of 100 ms. Individual images have a 
S/N ratio of 250 in the continuum. For each position of the slit on 
the sun, 360 images were recorded and summed to improve the 
S/N ratio to about 5000. The curvature of lines was corrected by 
a parabolic fit derived from the strong Fe i 630.15 nm line close 
to the [OI] 630 nm line with flat-held data. Then, a summation 
was done in the spatial direction of the slit to get high-precision 
spectra for various yu-values with a final S/N ratio of 100000 in 
the continuum. 

Eight positions on the solar disc are observed, at = 
{1.00,0.96,0.88,0.72,0.66,0.48,0.35,0.25}. 


4. Data analysis 

In our previous analysis of the solar oxygen abundance based 
on high-resolution, high S/N, central intensity and flux spec¬ 
tra (Paper I), we found that the oxygen abundance A(0) 
from the 630 nm line is 0.1 dex lower than the value obtained 
from the 636 nm line. It is known that the 630 nm line is 
blende d with a Nil line, wh ose \oggf has been well deter¬ 
mined (l.lohansson et al.ll2003l) . On the other hand the 636 nm 
line, which is weaker than the other [OI] line, lies on the red wing 
of a Ca I autoionisation line, and is blended with CN molecular 
lines, whose \og gf are more uncertain. From these considera¬ 
tions one should give preference to the A(0) determination from 
the 630 nm line, and for this reason we concentrate on this line 
in the present work. 

The two forbidden lines are very similar from the point of 
view of atomic data, and are formed in the same atmospheric 
layers. Based on the 3D model, we computed the contribution 
function for the equivalent width (EW) of the [OI] 630 nm line 
at the various p-values of the Kitt Peak intensity spectra with the 
spectral synthesis code Linfor3D. These are shown in Pig.[Tl As 
yU decreases, the contribution function becomes broader and the 
peak moves to higher photospheric layers. Eor // = 1.0 the max¬ 
imum value of the contribution function is at log Teao nm = -0.6 
and has a full width at half maximum (EWHM) of 1.2 in log r; 
for yu = 0.25 the maximum is at -1.0 and the EWHM 1.5 in log t. 
The spectra at different limb angles sample slightly different 
thermodynamic conditions at which the 630 nm blend is formed. 
In principle, different spectral lines probe different atmospheric 
layers, and show a different centre-to-limb variation. This can 
help to disentangle the contribution of oxygen and nickel to the 
the total equivalent width of the blend. In Eig.|2] the contribu¬ 
tion function of the EW is shown for the Ni contribution. The 
similarity of the yu-dependence of the contribution functions of 
oxygen (Eig.[TJ and nickel (Eig.[2]i already indicates that it is still 
a difficult task. 

The feature we investigate is a O-Ni blend. The observed 
spectra at disc-centre agree in that the EW of the blend is of 
the order of 0.47 pm. It is clear that the derived abundances of 
oxygen and nickel are related: any increase in A(0) implies a 
decrease in A(Ni) and vice-versa. A high oxygen abundance has 
to be associated with a low Ni abundance (a high A(0) implies 
a high contribution to the EW of the blend due to oxygen so that 
the fraction of EW left for Ni is small, implying a low A(Ni)) 
and a low A(0) to a high A(Ni). The situation is summarised in 
Fig .|3] where we present the relative abundances of oxygen and 
nickel for a given range in EW based on 3D analysis. 


p = 0.10,0.12,0.17,0.25,0.35,0.48,0.66,0.87,1.00 



Fig. 1. The contribution function of the EW of the pure [OI] 
630 nm line for various //-angles as a function of monochromatic 
optical depth in the continuum. The //-values are given at the 
top of the panel, the curves are plotted alternating as solid and 
dashed lines. The ones peaking at higher optical depth corre¬ 
spond to higher //-values. 


p = 0.10,0.12,0.17,0.25,0.35,0.48,0.66,0.87,1.00 



log X03O 

Fig. 2. Like Eig.[T]but for the blending Ni i line at 630 nm. 

4. 1. Line profile fitting 

We determine the oxygen abundance of the [OI] line at 630 nm 
from line profile fitting. In this way we take into account 
the line strength and the line shape at the same time. As in 
iBonifacio & Caff^ (l2003l) . we fitted the line profile interpolat¬ 
ing in a grid of synthetic spectra. To use the complete informa¬ 
tion, we simultaneously fitted all the limb angles of the solar 
intensity profiles, discarding the spectra closest to the limb. The 
same method is also us ed to analyse the O i triplet lines and is 
extensively explained injSteflfen,et_aD^2015). The previous anal¬ 
ysis of this [OI] line bv iPereira et al.l (l2009b I focused only on the 
centre-disc spectrum for the line profile fitting procedure. 

To fit the observed spectra with the grid of 3D spectra, we 
must normalise the wings of the [OI] line to unity because the 
3D synthetic spectra contain only the lines of oxygen and nickel. 
No molecular or other line contributions that would prevent nor¬ 
malised spectra from going to unity are considered in the spec¬ 
trum synthesis. To try to be as objective as possible, and to avoid 
a subjective normalisation of the spectra, we select two ranges 
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A(0) 

Fig. 3. The plot shows the oxygen and nickel abundances derived 
in the past years. The grey bar is the presently discussed range of 
the solar Ni abundance. Red symbols mean disc-integrated spec¬ 
trum, green symbols disc centre. Stars indicate theoretical mod¬ 
els, squares the HM model. The round dot of CentenoOS: in the 
work a sun-spot spectrum is used. The dashed lines are lines of 
constant EW for disc centre, solid for disc-integra ted light. This 
figure is an updated version of the one shown by iLudwig et alJ 
(l2010l) . The big green-blue star represents the result of this inves¬ 
tigation. The horizontal dotted line a t A(Ni)=6.23 represents val¬ 
ues from iGrevesse & Sauval (Il998h: the horizontal d ot-dashed 
line at A(Ni)=6.17 the value from iscott et al. I (12009 ): the hori- 
zontal dot-dashed line at A(Ni)=6.20 the value from IScott et alJ 
(1201 .^h . 



6300.2 6300.3 6300.4 

Wavelength (A) 


Fig. 4. VTT observations with a resolving power of 700000 
(black circles) compared to the 3D synthetic profile (solid red) 
at p = 0.4 (top) and at disc-centre (bottom). Below each profile 
the residual (xlO, displaced by h- 0.94) is shown. 


that are not contaminated by lines on the red and on the blue side 
of the [OI] line, and use these two ranges to pseudo-normalise 
the spectra at all //-angles. Using this procedure, there are be¬ 
tween 10 and 18 free fitting parameters (A(0), A(Ni), and wave¬ 
length shift and continuum level for each of the spectra, observed 
at four to eight //-values). We find that the fits are not always very 
robust. The minimum of the is not narrow and in the region 
of the best A(0) and A(Ni) values there are several local min¬ 
ima. By changing the input parameters, the values we derive for 
A(0) and A(Ni) are not always the same but fall in a limited 
range spanning ^ 0.02 dex. We restricted the fit to reasonable 
values of oxygen abundance (8.6 < A(0) < 8.85) and nickel 
abundance (5.85 < A(Ni) < 6.35). 

Another problem that needs to be mentioned is that we did 
not apply any instrumental broadening to the synthetic profile 
(which would correspond to about 0.5km/s for the VTT data 
and 1.5 km/s for Hinode spectrum) because it is already as broad 
as or slightly broader than the observed one even without any ex¬ 
tra broadening. The problem is visible in Fig.|4]where two VTT 
profiles, observed at centre-disc and at // = 0.4, respectively, are 
compared to a synthetic profile with the oxygen and nickel abun¬ 
dances of A(0)=8.72 and A(Ni)=6.17. From the figure it is clear 
that the agreement is worse towards the limb. 

The five sets of observations provide A(0) and A(Ni) in 
limited ranges. The results are summarised in Table|2l When 
we average the results we obtain A(0) = 8.735 + 0.017 and 
A (Ni) = 6.096 + 0.049. If we restrict our average to the data 
of Hinode, VTT, and SST, which we think are of higher quality, 
we obtain: A(0) = 8.725 + 0.012 and A(Ni) = 6.122 + 0.034, 
the uncertainty being the observed-to-observed scatter. This lat¬ 
ter value for the Ni abundance is in agreement with the value of 


Table 2. Summary of the results of the line-profile fitting. 


Spectrum 

A(0) 

A(Ni) 

Hinode 

8.71 

6.16 

SST 

8.73 

6.13 

WCL 

8.74 

6.10 

VTT 

8.74 

6.08 

THEMIS 

8.76 

6.02 


A(Ni) = 6.17 + 0.02 fromIScott et al. I (12009) and in reasonable 


agreement with A(Ni) = 6.20 + 0.04 from IScott et al.l (1201 5h . 
w ho used an improv ed 3D solar model atmosphere with respect 
tolScott et al. I(l2009t). Our nickel abundance is almost within one 


cr of the Scott et akT (l2009h and of the lScott et al.l (1201 5l) results. 


4.1.1. Hinode data 

We investigated the possibility of scattered light, but it does not 
seem to be a problem for the Hinode observations. The compar¬ 
ison of the Hinode centre-disc spectrum with the Neckel centre- 
disc data gives a perfect agreement (see Fig.|5]l. 

We fitted the five observed spectra corresponding to //-values 
of 0.99, 0.95, 0.86, 0.70, and 0.38 for two sets of Hinode data. 
The best fit of the five //-values gave A(0)=8.71 and A(Ni)=6.16 
(see Fig. |6li. 
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Fig. 5. Comparison of the Hinode centre-disc observation 
(dashed red) with the Neckel centre-disc spectrum (solid black); 
the agreement is very good. The Neckel spectrum has been 
broadened by 2 km/s for comparison with the Hinode observa¬ 
tion. In the lower panel the residual (xlOO) is shown. 


4.1.2. SSTdata 

The best fit of the four p-values gave A(0)=8.73 and 
A(Ni)=6.13 (see Fig.|2l). The fit looks good and the result is in 
agreement with the values derived from the other observed spec¬ 
tra. The spectra show no signs of scattered light. 

4.1.3. WCL data 

We compare d the WCL obse r ved sp ectr a at ^ = 1.0 to t he so- 
lar atlases of iDelbouille et all (119731) and iNeckel & Lab^ (Il984t) 
and discovered that our WCL disc-centre spectrum has a smaller 
EW (by roughly 4%) than in both high-resolution solar atlases. 
This could be due to scattered light in the spectrograph, solar 
activity, or some other unknown effect. We corrected the data by 
this factor, assuming that the scattered light affects the spectra at 
all p-angles in the same way. We restricted the line profile fitting 
to the six spectra with fi > 0.2. Using this procedure, we fitted 
the observed spectra with 14 free parameters (A(0), A(Ni), and 
wavelength shift and continuum level for each of the six spectra) 
and obtained A(0)=8.74 and A(Ni)=6.10 (see Fig.|3). 

4.1.4. VTT data 

We find no evidence for a contamination by scattered light in 
these spectra. 

We note that we are not able to fit the high-quality VTT spec¬ 
tra because our synthetic profiles are too broad. We therefore 
decided to broaden the observed profiles and we were able to 
obtain a very good fit (see Fig.0for an example). The required 
extra broadening is //-dependent: 1.3, 1.7, 1.8, and 1.9km s“* 
for // positions 1.0, 0.8, 0.6, and 0.4, respectively. In this way we 
derive A(0)=8.74 and A(Ni)=6.08. With five //-values, the most 
extreme limb component, // = 0.3, required an extra broadening 
of 1.9 kms“' . In this case we derive A(0)=8.72 and A(Ni)=6.12 
(see Fig. 12). The oxygen and nickel abundances we derive in the 
two cases are very similar, but we prefer to use the case with four 
//-values because the agreement with the observed data is better. 



6300. 1 6300. 2 6300. 3 6300. 'I 

uauelength (A) 


Fig. 6. Fitted 3D line profiles for the five //-values (red) com¬ 
pared to the Hinode observed spectra (black). From the simul¬ 
taneous fit we derive A(0)=8.71 and A(Ni)=6.16. Below each 
profile the residual (xlO, displaced by -1-0.94) is shown. 


4.1.5. THEMIS data 

As we did in the previous cases, we compared the observed 
line profile of the disc-centre spectrum to the solar atlases of 



Fig. 7. Fitted 3D line profiles for the four //-values (red) com¬ 
pared to the SST observed spectra (black). From the simultane¬ 
ous fit we derive A(0)=8.73 and A(Ni)=6.13. Below each profile 
the residual (xlO, displaced by -(-0.94) is shown. 
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6300.1 6300. 2 6300. 3 6300 .4 6300.1 6300. 2 6300.3 6300 .4 

Uavelength (A) Uauelength (A) 


Fig. 8. Fitted 3D line profiles for the six //-values (red) com¬ 
pared to the WCL observed spectra (black). From the fit we de¬ 
rive A(0)=8.74 and A(Ni)=6.10. Below each profile the residual 
(x5, displaced by -1-0.94) is shown. 

Table 3. 1000 Monte Carlo simulations on synthetic spectra. 


S/N 

A(0) 

A(Ni) 

300 

8.66 ± 0.04 

6.23 ± 0.08 

400 

8.66 ± 0.03 

6.24 ± 0.07 

800 

8.67 ± 0.02 

6.24 ± 0.04 


iDelbouille et al.l (Il973h and iNeckel & Labsl (Il984l) . From this 
comparison we estimate a scattered light fraction of about 4%. 
The quality of the THEMIS spectra is inferior to that of the other 
sets of observations, so that the continuum is difficult to define. 
This could be due to a not entirely successful flat-held correction 
of the data. From the fit shown in Fig.[T0]we obtain A(0)=8.76 
and A(Ni)=6.02, but the fit is not as good as for the other obser¬ 
vations. 

4.1.6. Monte Carlo tests 

We investigated the role of the S/N ratio on the uncertainty in 
the abundance determination. We performed two different Monte 
Carlo simulations. In the first, we injected Poisson noise in the 
synthetic spectra with parameters A(0)=8.66 and A(Ni)=6.25. 
We considered the five //-values of the Hinode data, 1.00, 0.95, 
0.86, 0.70, and 0.38. The results are presented in Tabled The 
uncertainties in the Table are fit-to-fit scatter for three cases of 
noise injection of signal-to-noise of 800, 400, and 300. We can 
consider these uncertainties as the statistical error in the abun¬ 
dance determinations related to Poisson noise. Owing to the 
high-quality spectra we work with, we suggest an uncertainty 
of 0.02 dex and 0.04 dex for the derived abundance of oxygen 
and nickel, respectively. 



G300.1 6300.2 6300.3 6300.4 

uave length (A) 


Fig. 9. Fitted 3D line profiles for five //-values (red) compared 
to the VTT observed spectra (black). From this fit we derive 
A(0)=8.72 and A(Ni)=6.12. Below each profile the residual 
(xlO, displaced by -1-0.94) is shown. 


We then took the observed Hinode spectra, slightly degrad¬ 
ing their quality by injecting noise, but obtaining a plausible dif¬ 
ferent “realisation” of an observation in order to verify whether 
the result is consistent with the previous exercise of the Monte 
Carlo simulation over synthetic profiles. This Monte Carlo sim¬ 
ulation consists of 1000 events where Poisson noise is injected 
into the five Hinode observed spectra. The results are presented 
in Table|4] We see that the scatter within the realisations for oxy¬ 
gen is twice the value derived from synthetic spectra, which is 
expected because we start from an observed spectrum. For Ni the 
fit-to-fit scatter is about a factor of three larger than in the case 
of the synthetic spectra. Nickel seems to be more difficult to de¬ 
rive from observed spectra because the line is weaker, or perhaps 
because there are unce rtainties in the wavele ngths of the Ni iso¬ 
topic components (see iRosberg et al.lll993l) . The results of the 
Monte Carlo tests give an indication that the systematics domi¬ 
nate the error budget. 


72 








































Caffau et al.: Solar oxygen project III 



6300.1 6300.2 6300.3 6300.4 6300.1 6300.2 6300. 3 6300.4 

Ua^^elength (A) Uavelength (A) 


Fig. 10. Fitted 3D line profiles at the eight //-values (red) com¬ 
pared to the THEMIS observed spectra (black). From the fit we 
derive A(0)=8.76 and A(Ni)=6.02. Below each profile the resid¬ 
ual (x5, displaced by h- 0.94) is shown. 

Table 4. 1000 Monte Carlo simulations on the Hinode spectra. 


S/N 

A(0) 

A(Ni) 

300 

8.71+0.06 

6.13 + 0.14 

400 

8.71+0.06 

6.14 + 0.14 

800 

8.70 + 0.05 

6.15 + 0.13 


5. Discussion 

5. 1. The blending nickel line 

In the previous sections we have discussed the role of the nickel 
blend in terms of the Ni abundance. However, we believe it is 
unlikely that a re-analysis of the solar Ni abundance will result 
in a la rge downward revision. The recent result of IScott et al.l 
(l2015ll . who find A(Ni) = 6.20 + 0.04, corroborates this no¬ 
tion. For the meteoritic nickel abundance, iLoddersI (l2003h rec¬ 
ommends A(Ni)= 6.22 + 0.03, thus a nickel abundance as low as 
around 6.00 (THEMIS) would be in stark disagreement with the 
meteoritic nickel abundance. Likewise, such a revision cannot be 
due to errors in the /-value of the Ni i line, which i s known with 
an uncertainty of only 14% (I.Tohansson et al.ll2003h . i.e. 0.06 dex 
in logarithmic abundance. On the other hand, we note that the 
Ni line has a rather high excitation energy (4.266 eV) and that 
no calculations exist concerning the deviations from LTE for 


Ni. NLTE effects may produce a line weaker than p redicted by 
our LTE computations. However. IScott et al. I (l2009h obtained a 
rather homogeneous Ni abundance in LTE from all Nil lines 
used in their analysis, which would make strong departures from 
LTE for the blending Ni line rather surprising. Computations on 
a Fe I line with similar atomic properties provided by Lyudmila 
Mashonkina (priv. comm.) suggest that the NLTE correction of 
the Ni I line cannot be the solution to the problem. 

We consider here a Ni abundance of about 6.1 dex, which - 
taking into account the uncertainty in the o scillator strength of 
the Ni line - is perfectly consistent with the IScott et ^ (120091) 
result. We would like to emphasise that a low contribution of 
nickel to the 630 nm blend corresponding to A(Ni) ss 6.1 is 
also cle arly compatible with previo us findings. Combining the 
result of lAllende Prieto et al.l (1200 ll) . logg/eNi = 3.94, with the 
presently accepted value of the oscillator strength (logg/Ni = 
-2.11) results in a nicke l abunda nce of 6.05. We n ote that 
lAllende Prieto et ^ (1200 ll) and also IScott et al. I (l2009l) applied 
the same 3D model atmosphere in their analyses, but Scott et 
al. arrived at a sig nificantly higher nickel abund ance from the 
unblended Ni lines. iMelendez & Asplundl (l2008l) used the [OI] 
line at 557.34nm to obtain an oxygen abundance of 8.71 + 0.07, 
which is in reasonable agreement with our res ult for t he two 
other forbidden lines within the uncertainties. lAvres I (l2008l) 
analysed the 630 nm blend with a single snapshot from the 
same solar CO^BOLD model as used in this work. He obtains 
A(0) = 8.81+0.04 and A(Ni) = 6.08 (he states a ^ 30 % s maller 
Ni abundance than the value given bv iGrevesse & Sauvall d 19981) 
which is 6.23 ± 0.04) adopting a rather large equivalent width for 
the feature at disc-centre of 0.476 pm, and applying corrections 
to the model to reproduce the measured continuum intensity at 
the wavelength of the blend. While his result is certainly at the 
high end at the presently discussed oxygen abundances, again, 
the contribution to the blend by nickel is small. 

5.2. The role of the centre-to-limb variation 

One should bear in mind that the purpose of this work was to 
investigate the line formation in the solar photosphere at differ¬ 
ent limb-angles. In this respect the data we used are unique, but 
not all these spectra have the high quality of the traditional solar 
atlases. As we pointed out above, the equivalent widths (EWs) 
of the two WCL and THEMIS disc-centre spectra of the 630 nm 
line are smaller than the ones of the largely used disc- centre solar 
atlases (iDelbouille et al.l[l973l : lNeckel & Labslll984l) . However, 
we argue that our observations in their ensemble are surely suffi¬ 
cient to allow a differential study of the centre-to-limb variation 
of the forbidden oxygen feature. 

As pointed out in Section l4.1 .41 the synthetic profiles are too 
broad when compared to high-quality data of extremely high 
spectral resolution. The extra broadening in the synthetic spectra 
could indicate a velocity field in the hydrodynamical model that 
is too strong. This effect is most evident in the VTT observation 
owing to their high resolution, but the effect that the FWHM of 
the observed spectrum becomes increasingly smaller than in the 
synthesis towards the limb is visible in the SST (see Fig.|7]i and 
THEMIS data as well (see Fig.[l0]) 

6. Conclusions 

Our results can be summarised as follows: 

1. The oxygen abundance we derive from the 630 nm blend is 

A(0)=8.73 ± 0.02 ± 0.05 (with error estimates from the scat- 
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ter of the different observations and the Monte Carlo simula¬ 
tions). 

2. The above oxygen abundance from the 630 nm feature is in¬ 
consistent with a Ni abundance of 6.25 and the assumption 
of LTE: we advocate a lower contribution of the Ni line to the 
blend corresponding to a nickel abundance of 6.11 + 0.04, or 
logg/eNi =4.00. 

3. The synthetic line profile we obtain from the 3D model 
is broader than expected, and this is particularly evident 
in comparison with the VTT observations with the highest 
spectral resolving power of 700000. This effect, which is not 
seen for other lines, can be due to an overly large velocity 
field in the model atmosphere, although this effect may also 
indicate errors in the adopted atomic line wavelength differ¬ 
ence between the two components of the [OI]-i-Ni i blend. 

4. The observed spectra at the lowest limb angles (// < 0.2) are 
not reproduced satisfactorily by our computations. We spec¬ 
ulate that this problem is due to the solar model that is neither 
taking into account sphericity effects nor a chromosphere. 

5. On the other hand, the synthetic spectra reproduce most of 
the observed data for fi > 0.2 very well, even though the 
agreement between observed and synthetic spectra deterio¬ 
rates somewhat towards the limb. 

6 . With the low oxygen abundance we derive in this work from 
the [OI] line at 630 nm, the discrepancy of about 0.1 dex with 
the oxygen abundance derived from the [OI] 636 nm line in 
Paper I and discussed in Paper II is reduced to 0.06 dex due 
to a lower A(Ni), but remains unexplained. 
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